Explosive Synchronization Transitions in Scale-free Networks 
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The emergence of explosive collective phenomena has recently attracted much attention due to the discovery 
of an explosive percolation transition in complex networks. In this Letter, we demonstrate how an explosive 
transition shows up in the synchronization of complex heterogeneous networks by incorporating a microscopic 
correlation between the structural and the dynamical properties of the system. The characteristics of this ex- 
plosive transition are analytically studied in a star graph reproducing the results obtained in synthetic scale-free 
networks. Our findings represent the first abrupt synchronization transition in complex networks thus providing 
a deeper understanding of the microscopic roots of explosive critical phenomena. 



PACS numbers: 



l.20.-a. 



I.75.HC, 89.75.Kd 



o 



> 
m 

00 
O 



X 



Synchronization is one of the central phenomena represent- 
ing the emergence of collective behavior in natural and syn- 
thetic complex systems lU-lstl- Synchronization processes de- 
scribe the coherent dynamics of a large ensemble of intercon- 
nected autonomous dynamical units, such as neurons, fireflies 
or cardiac pacemakers. The seminal works of Watts and Stro- 
gatz Hut] pointed out the importance of the structure of inter- 
actions between units in the emergence of synchronization, 
which gave rise to the modern framework of complex net- 
works IS]. 

Since then, the phase transition towards synchronization 
has been widely studied by considering non-trivial networked 
interaction patterns ITJ. Recent results have shown that the 
topological features of such networks strongly influence both 
the value of the critical coupling, Ac, for the onset of synchro- 
nization H - 12 ] and the stability of the fully synchronized state 
113-3. The case of scale-free (SF) networks has deserved 
special attention as they are ubiquitously found to represent 
the backbone of many complex systems. However, the topo- 
logical properties of the underlying network do not appear to 
affect the order of the synchronization phase transition, whose 
second-order nature remains unaltered [jsl . 

More recently, the study of explosive phase transitions in 
complex networks has attracted a lot of attention since the 
discovery of an abrupt percolation transition in random [llTIl 



and SF networks II18L 11911 . However, several questions about 



the microscopic mechanisms responsible of such an explosive 
percolation transition and their possible existence in other dy- 
namical contexts remain open. In this line, we conjecture that 
such dynamical abrupt changes occur when both, the local 
heterogeneous structure of networks and the dynamics on top 
of it, are positively correlated. 

In this Letter, we prove our conjecture in the context of the 
synchronization of Kuramoto oscillators. We show that an 
explosive synchronization transition emerges in SF networks 
when the natural frequency of the dynamical units are posi- 
tively correlated with the degree of the units. Furthermore, 
we analytically study this first-order transition in a star graph 



and show that the combination of heterogeneity and the above 
correlation between structural and dynamical features are at 
the core of the explosive synchronization transition. 

Let us consider an unweighted and undirected network of 
N coupled phase-oscillators. The phase of each oscillator, 
denoted by Oi{t) (i = 1, ■ • ■ , N), evolves in time according to 
the Kuramoto model 120(1 : 
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withi = 1,...,^ (1) 



where uji stands for the natural frequency of oscillator i. The 
connections among oscillators are encoded in the adjacency 



matrix of the network. A, so that Ai 



1 when oscillators 



i and j are connected while Aij — otherwise. Finally, the 
parameter A accounts for the strength of the coupling among 
interconnected nodes. 

The original Kuramoto model assumed that the oscillators 
were connected all-to-all, i.e. Aij = 1 Vi ^ j. In this setting, 
a synchronized state, i.e. a state in which 9i{t) = Oj{t) 
and Vt, shows up when the strength of the coupling A is larger 
than a critical value To monitor such synchronization 

transition as A grows, the following complex order parameter, 
which quantifies the degree of synchronization among the N 
oscillators, is used 112 311 : 
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The modulus of the above order parameter, r{t) e [0, 1], 
measures the coherence of the collective motion, reaching the 
value r = 1 when the system is fully synchronized, while 
r = for the incoherent solution. On the other hand, the value 
of ^'(t) accounts of the average phase of the collective dynam- 
ics of the system. Typically, the average (over long enough 
times) value of r as a function of the coupling strength A dis- 
plays a second-order phase transition from r = to r = 1 
with a critical coupling A^ = 2/{TTg{Ld = 0)), where giuj) 



2 





. ■■' 






Forward — »— f 




Backward — • — 4 




■ 


(a) 





0.8 
0.6 



f Forward 
f Backward 



^ 



0.2 0.4 0.6 0.8 1 1.2 0.4 0.6 0.8 1 1.2 1.4 



0.8 1 1.2 1.4 1.6 1.8 0.8 1 1.2 14 1.6 1.8 

FTG. 1: (color online) Synchronization diagrams r(A) for different 
networks constructed using the interpolation model introduced in 
(2^ . The a values in each panel are (a) a = 1 (ER), (b) a = 0.6, (c) 
a — 0.2 and (d) a = (BA). The four panels show both Forward 
and Backward continuations in A using increments of 5\ = 0.02. 
The size of the networks is A*" = 10'' and the average degree is 
{k) = 6. 



is the distribution of the natural frequencies, {wi}, and it is 
assumed to be unimodal and even 12311 . 

Here we will focus on the influence of the dynamical and 
topological characteristics at the local level (the nodes of the 
network and their interactions) in the emergence of global 
synchronization. In particular, we will identify the internal 
frequency of each node i directly with its degree ki, so that 
uji = ki in Eqs. ([U. Note that this prescription automatically 
sets that the distribution of frequencies g{uj) = P{k) but not 
vice versa 1I24I1 . 

To study the effects of the correlation between dynamical 
and structural attributes, we simulate the Kuramoto model 
on top of a family of networks generated according to 1125 1. 
This model allows to construct networks with the same av- 
erage connectivity, (fc), interpolating from Erdos-Renyi (ER) 
graphs to Barabasi-Albert (BA) SF networks by tuning a sin- 
gle parameter a. The growth of the networks assumes that a 
newly added node either attaches randomly with probability a 
or preferentially to those nodes with large degree with proba- 
bility (1 — a). In this way, a = 1 gives rise to ER graphs with 
a Poissonian degree distribution whereas for a = the result- 
ing networks are SF with P{k) k~^. Intermediate values 
a G (0, 1) tune the heterogeneity of the network, which in- 
creases when going from a = 1 to a = 0. In the four panels 
of Fig.[Tl we report the synchronization diagrams of four net- 
work topologies constructed using this model. The limiting 
cases of ER and BA networks correspond to panel [T^ and[TJl, 
respectively. The size of these networks are = 10'^ while 
the average connectivity is set to {k) = 6. 

For each panel in Fig. [T] we have computed two synchro- 
nization diagrams, r(A), labeled as Forward and Backward 
continuations. The former diagram is computed by increas- 



ing progressively the value of A and computing the stationary 
value of the order parameter r for Aq, Aq + S\,..., Xq + nSX. 
Alternatively, the backward continuation is performed by de- 
creasing the values of A from Xq + n6X to Aq. The panels [T^, 
[T]5 and[T]; show a typical second-order transition with a perfect 
match between the backward and forward synchronization di- 
agrams. Importantly, the onset of synchronization is delayed 
as the heterogeneity of the underlying graph (and thus that of 
the frequency distribution g{uj)) increases. 

The most striking result is however observed for the BA 
network (panel [TJl) in which a sharp, first-order synchroniza- 
tion transition appears. In the case of the forward continuation 
diagram the order parameter remains r ~ until the onset of 
synchronization in which r jumps suddenly to r ~ 1 pointing 
out that almost all the network has reached the synchronous 
motion. Moreover, the diagram corresponding to the back- 
ward continuation also shows a sharp transition from the fully 
synchronized state to the incoherent one. The two sharp tran- 
sitions takes place at different values of r so that the whole 
synchronization diagram displays a strong hysteresis. 

To analyze deeply the change of the order of the synchro- 
nization transition, we have monitored the evolution of the dy- 
namics for every node by computing their effective frequency 
along the forward continuation, see Fig.|2]. The effective fre- 
quency of a node i is defined as 
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0^{T) dr , 



(3) 



with T ^ 1. We have also computed the evolution of wf^ 
within a degree class k, {uj)k, averaging over nodes having 
identical degree k: 
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where Nk — NP{k) is the number of nodes with degree k 
in the network. From the panels in Fig.|2]we observe that the 
individual frequencies and the different curves {uj)k{X) con- 
verge progressively to the average frequency of the system 
= (fc) = 6 until full synchronization is achieved. Panel |2^ 
(ER graph) shows that the convergence to is first achieved 
by those nodes with large degree while the small fc-classes 
achieve full synchronization later on. As the heterogeneity of 
the network increases (see a = 0.6 and a = 0.2 in panels 
I2J3 and 12};, respectively) the differences in the convergence of 
the fc-classes decrease. Finally, for the BA network (Fig. |2jJ), 
we observe that nodes (and thus the different fc-classes) re- 
tain their natural frequencies until they become almost fully 
locked, which signal the abrupt synchronization observed in 
Fig. [TJl. Thus, the first-order transition of the BA network 
corresponds to a process in which no microscopic signals of 
synchronization are observed until the critical coupling Ac is 
reached. 

To further explore the correspondence of the explosive syn- 
chronization transition with the SF nature of the underlying 
graph, in Fig. [3]a we show the synchronization diagrams for 
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FIG. 2: (color online) The panels show the evolution of the effective 
frequencies of the nodes along the (forward) continuation in the four 
model networks of Fig. [T] The colored dots account for single-node 
values (colors stand for their respective degree) while the solid lines 
show the evolution of the average value of the effective frequencies 
of nodes having the same degree. 



different uncorrelated SF graphs with different degree distri- 
bution' exponents. These graphs have been constructed using 
the configurational model ll26ll by imposing a degree distri- 
bution P{k) - fc-T with 7 = 2.4, 2.7, 3.0 and 3.3. The 
synchronization diagrams are obtained by forward continua- 
tion (as described above) starting at A = 1 and performing 
adiabatic increments of b\ = 0.02. Again, for each value of 
A the Kuramoto dynamics is run until the value of r reaches 
its stationary state. From the figure it is clear that a first-order 
synchronization transition appears for all the reported values 
of 7 pointing out the ubiquity of the explosive synchroniza- 
tion transition in SF networks. Moreover, the onset of syn- 
chronization, Ac, is delayed as 7 decreases, i.e. when the het- 
erogeneity of the graph increases. 

Up to now, we have shown that the explosive synchroniza- 
tion transition appears in SF when the natural frequencies of 
the nodes are correlated with their degrees. To show that this 
correlation is the responsible of such explosive transition, in 
Fig. [Sji we show the synchronization diagram for the same 
SF networks used in Fig. [3^, but when the correlation be- 
tween dynamics and structure is broken in such a way that the 



same distribution for the internal frequencies, giiS) = uj~'^ is 
kept. To this end, we made a random assignment of frequen- 
cies to nodes according to 17(0;). The plots reveal that now 
all the transitions turn to be of second-order, thus recovering 
the usual picture of synchronization phenomena in complex 
networks. Therefore, the first-order transition arises due to 
the positive correlation between natural frequencies and the 
degrees of the nodes in SF networks ll27ll . 

All the simulations results presented corroborate our con- 
jecture about the explosive percolation transition in SF net- 
works. To get analytical insights, we reduce the problem stud- 
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FIG. 3: (color online) Panel (a) shows the synchronization diagrams 
r(A) for several SF networks constructed via the configurational 
model. All the networks have a degree-distribution P{k) ~ k ' 
with 7 = 2.4, 2.7, 3.0, and 3.3 while iV = 10^ The steps of 
the continuation are set to 5 A = 0.02. In panel (b) we show the 
synchronization diagrams of the same SF networks without the local 
correlation between degrees and natural frequencies, i.e. uji 7^ ki, 
while the distribution of natural frequencies is still P((jj) ~ a; 



ied to the analysis of the star configuration, a special structure 
that grasp the main property of SF networks, namely the role 
of hubs. Therefore, we explore the synchronization transition 
of such a configuration and show that it is indeed explosive 
when the correlation LOi = ki holds. A star graph (as shown in 
the inset of Fig. |4^) is composed by a central node (the hub) 
and K peripheral nodes (or leaves). Each of the peripheral 
nodes connects solely to the hub. Thus, the connectivity of 
the leaves is fc^ = 1 (z = 1, ...,K) while that of the hub is 
kh = K. Let us suppose that the hub has a frequency LOh 
while all the leaves beat at the same frequency w. 

First we set a reference frame rotating with the average 
phase of the system, '^{t) = *(0) + fit, being the av- 
erage frequency of the oscillators in the star, = {Kuj + 
ujfi)/{K + 1). In the following we set ^'(0) = without loss 
of generality so that the transformed variables are defined as 
<j)h = Qh-^t for the hub and 4>j = 6j - fit (with j = l,...,K) 
for the leaves. Thus, the equations of motion for the hub and 
the leaves read; 

K 

(fh = (w/i - f^) + A ^ sin(0j ~(j)h), (5) 

(pj = {lj — il) + Xsm{(f)h — (fij), with j = l...K.(6) 

In this rotating frame the motion of the hub, Eq. (|5]l, can be 
expressed as: 



<^h = (w/i -n) + \{K + l)rsvcL{(l>h) 



(7) 



note that in this new frame it is easy to identify that the dynam- 
ics of the hub is governed by its new inherent frequency and 
the superposition of a set of identical signals from the leaves. 
Now, imposing that the phase of the hub is locked, (j)h = 0, 
we obtain: 



A(/i + l)r 



(8) 
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FIG. 4: (color online) We show the synchronization diagrams for the 
star graph [see the inset in plot (a)]. In (a) we show the (forward 
and backward) continuation diagrams for the case K = 10 while (b) 
shows the forward continuation diagrams for different star graphs of 
different sizes, corresponding to K = 20, 30, 40, 50, 60 and 70. 



of the distribution of natural frequencies. Our findings pro- 
vide with an explosive phase transition of an important macro- 
scopic phenomena, synchronization, in a widely studied dy- 
namical framework, the Kuramoto model, thus shedding light 
to the microscopic roots behind these phenomena and paving 
the way to their study in other dynamical contexts. 
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Now we consider the equations for the leaves, Eq. (|6]l, and 
evaluate the expression for cos (pj in the locked regime, (pj = 
0. After some algebra, we obtain the following expression; 



{n ~ uj) sin(t>h ±J[^- sin^ [A^ - {n - uj)^] 

cos = 

A 

(9) 

The above expression is valid only when {ft — u) < A. From 
this inequality we obtain the value of the coupling A for which 
the phase-locking is lost, i.e. the critical coupling Ac = il — ui. 
In our case, we have Uh = K, lj = 1 and il = 2K/ {K + 1) 
so that we obtain a critical coupling Ac = {K — 1)/ {K + 1). 
On the other hand, we can derive the value Tc of the order 
parameter at the critical point by using Eq. (|8} and Eq. (|9} to 
compute r ~ (cos(0)) at Ac: 



cos{(l)h) + Kcos{(l)j) 



K + 1 



K 



{K + 1) 



(10) 



Therefore, as Tc > 0, when the synchronization is lost there 
is a gap in the synchronization diagram pointing out the ex- 
istence of a first-order synchronization transition. Moreover, 
as K increases both Ac and Vc tend to 1 thus confirming the 
first-order nature of the transition in the thermodynamic limit, 
K ^ oo. As shown in Fig.|4^ for the case if = 10, the the- 
oretical values of Ac and Tc, are in perfect agreement with re- 
sults from numerical simulations. Finally, as shown in Fig.l^', 
the stability of the unlocked phase regime, r ~ 0, increases 
with K so that we can reach larger values of A by continuing 
(forward) the solution with r ~ 0. As a result, this robustness 
of the incoherent solution leads to a larger hysteresis cycle of 
the synchronization diagram as K increases. 

Summing up, we have shown that an explosive synchro- 
nization transition occurs in SF networks when there is a pos- 
itive correlation between the structural (the degrees) and dy- 
namical (natural frequencies) properties of the nodes. This 
constitutes the first example of an explosive synchronization 
transition in complex networks. Moreover, we have shown 
that the emergence of such transition is intrinsically due to the 
interplay between the local structure and the internal dynam- 
ics of nodes rather than being caused by any particular form 
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